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We show that semiconductor nanowires coupled to an s-wave superconductor provide a playground to study 
effects of interactions between different topological superconducting phases supporting Majorana zero-energy 
modes. We consider quasi-one dimensional system where the topological phases emerge from different trans- 
verse subbands in the nanowire. In a certain parameter space, we show that there is a multicritical point in the 
phase diagram where the low-energy theory is equivalent to the one describing two coupled Majorana chains. 
We study effect of interactions as well as symmetry-breaking perturbations on the topological phase diagram in 
the vicinity of this multicritical point. Our results shed light on the stability of the topological phase around the 
multicritical point and have important implications for the experiments on Majorana nanowires. 



I. INTRODUCTION. 

The possibility of realizing Majorana fermions, elusive 
particles that are their own anti-particles, in semiconductor 
nanowires coupled with an s-wave superconductor has at- 
tracted a lot of attention recently ^ In addition to the in- 
trinsic motivation of finding Majorana particles in nature^, 
the solid-state Majoranas have additional property of funda- 
mental physics interest: Majorana zero-energy modes emerg- 
ing in topological superconductors obey non-Abelian braid- 
ing statistics^ ^ which can be exploited for quantum compu- 
tation purposes^. The prediction of the emergence of Majo- 
rana fermions in semiconductor/superconductor heterostruc- 
tures^ has led to much activity aimed at detecting these ex- 
otic particleJ^^EU as well as exploiting them for topological 
quantum computation 

There is no doubt that Majorana fermions can be reali- 
azed in suitable mean-field models describing realistic physi- 
cal systems. The existence of Majorana zero-energy modes in 
these system can be shown theoretically b y explicitl y solving 
the corresponding quadratic Hamiltoniang^ECSIIlEol or invok- 
ing topological invariants developed for noninteracting sys- 
tems^^. The situation is much more complicated, however, 
once interactions are included, and there are examples where 
interactions lead to the breakdown of the classification devel- 
oped for noninteracting systems^^. In this paper we study 
effect of interactions on the topological phase diagram us- 
ing a realistic model which describes multiband nanowires 
proximity-coupled to an s-wave superconductor. The effective 
theory naturally emerging in multiband nanowires is equiva- 
lent to the model of two coupled Majorana chains. Rather 
than being spatially dependent, the coupling between Majo- 
rana "chains" is controlled by external parameters such as 
magnetic field and chemical potential, and can be tuned in a 
given sample. Within this model, we characterize the effect of 
interparticle interactions as well as various other perturbations 
on the topological phase diagram. We find that interactions do 
not change the phase diagram at the qualitative level but lead 
to a non-trivial renormalization of the phase boundary. 

The paper is organized as follows. In Sec. |ll| we first in- 
troduce theoretical model for multi-band Majorana nanowires 
and discuss topological phase diagram assuming no mixing 



between different subbands. We identify multicritical points 
in the phase diagram and derive an effective model describ- 



ing the multi-critical points. In Sec. Ill we consider all band- 
mixing terms allowed by the symmetry of the model and an- 
alyze their effect on the phase diagram, in particular around 
multicritical points. In Sec.|lV]we add interparticle interac- 
tion and calculate the topological phase diagram of the system 
using bosonization and real-space RG techniques. Finally, we 
conclude in Sec. [V| with the summary of main results. 



II. THEORETICAL MODEL. 

The system we consider here consists of a semiconductor 
quantum well with dimensions <^ Ly <^ in contact 
with an s-wave superconductor, see Fig[T^. We assume that 
the confinement along the 2; -axis is very strong so that only 
the lowest subband with respect to the z-axis eigenstates is 
occupied, whereas the confinement along the ?/-axis is much 
weaker and only a few ^/-subbands can be populated. Assum- 
ing that electrostatic gradient due to the applied gate poten- 
tial is along 2:-direction, the single-particle Hamiltonian takes 
the usual form for the 2D semiconductor with the spin-orbit 
Rashba interaction (h = 1): 
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where m*, a and ji are the effective mass, the strength of spin- 
orbit interaction and chemical potential, respectively. The lat- 
ter can be controlled using the gate electrodes^^ . The last 
term in Eq. ^ describes the Zeeman term due to an ap- 
plied external magnetic field aligned along the x-axis, Vx = 
9sm1^bBx/2. One can notice that Hamiltonian ^ is in- 
variant under transformation i/jx i%lj-\{—x) with A be- 
ing an electron spin. This transformation corresponds to a 
TT -rotation around x-axis combined with the inversion opera- 
tion P • /(x) = /(-x): 



UHU-^ = H where U = ia^P. 



(3) 



The proximity-induced superconductivity due to the ex- 
change of electrons between the semiconductor and supercon- 
ductor leads to the emergence of anomalous correlations in the 
semiconductor which are captured by the following Hamilto- 
nian: 



h.a] (4) 



We choose the gauge where the induced pair potential Aq is 
real. If Aq is x-independent which is assumed throughout the 
paper, the Hamiltonian ^ also respects symmetry ([3]). Thus, 
all terms that are derived from the total Hamiltonian Ht = 
H + Hsc should respect the symmetry ([3]). 

We now construct effective two-band model for the semi- 
conductor nanowire which captures the physics we are inter- 
ested in. First of all, the Hamiltonian ^ is separable inx — y 
coordinates and one can write the field operator as 



i^x{x,y)-- 
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where ax^p^^riy is electron annihilation operator in a state riy 
having momentum px and spin A. Assuming that the con- 
finement energy along ^/-direction E^h is larger than all the 
relevant energy scales of the Hamiltonian ([2]), one can project 
the Hamiltonian Ht to the lowest two subbands. Within this 
approximation, the Hamiltonian of the system reads (h = 1): 



dx 



-L 
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J W 
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dx 
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= J [Aoctc; + Aoc/t^^i + h.c] . (7) 



Here cx and dx represent fermion annihilation operators of the 
first and second subbands having spin A. The spin-orbit band 

('12) 

mixing Hamiltonian H^o originates from taking the expec- 
tation value py operator between different band eigenstates. 

The energy = dyf^ sin(^)a, sin(^) = ^. 
One can notice that a constant magnetic field does not 
lead to inter-subband mixing because the corresponding ma- 
trix elements vanish identically due to the orthogonality of 
the wavef unctions corresponding to different subbands. Thus, 
the two-band approximation is valid even for a large Zeeman 
splitting 14 - E^i,. 

It is implicitly assumed in Eqs. ([6])-(|7]) that the length of the 
wire 2Lx = 2L is much longer than the effective supercon- 
ducting coherence length ^ in the semiconductor. For simplic- 
ity, we take the induced superconducting pair potential to be 
the same for the first and second transverse bands. 




FIG. 1 . (Color online) a) Semiconductor quantum well with dimen- 
sions Lz <^ Ly <^ Lx in contact with an s-wave superconductor, b) 
Single-particle energy spectrum showing lowest two transverse sub- 
bands. The combination of the Rashba spin-orbit coupling and mag- 
netic field results in splitting of the spin degeneracy in each subband. 
c) Phase diagram for two-band topological superconductor assum- 
ing fermion parity in each subband is conserved. At the multicritical 
point the system has enhanced Z2 Z2 symmetry and is equivalent 
to the model of two Majorana chains. Changing of the mass terms 
Ma corresponds to a topological phase transition and allows one to 
map out the topological phase diagram. Here tilde denotes rescaled 
energy E = E/rrfa^ . d) Single-particle energy spectrum at the 
multicritical point where two bands belonging to different subbands 
touch, and the band topology changes in a non-trivial way. 



As shown in Ref . [24| the system described by the Hamilto- 
nian Hq = Hsu + realizes a non-trivial topological SC 
state in a suitable parameter regime. In the weak coupling 
limit Aq ^ 0, the topological phase emerges when there is an 
odd number of Fermi surfaces. In this case, the Hamiltonian 
of the system is adiabatically connected with the one of a spin- 
less p-wave superconductor which is known to host Majorana 
zero-energy modes at the ends of the nanowire. As shown 
in Fig. [TJ), such a situation is realized when \Vx\ > and 
I I > — ^sb I which corresponds to the topological phases 
originating from the first (Uy = 1) and second (Uy = 2) sub- 
bands, respectively. Looking at the single-particle band struc- 
ture as a function of Vx, one can notice that there is a special 
point Vx ~ Esh 1 2 where two bands belonging to = 1 and 
Uy = 2 subbands touch, see Fig.[T|: and d. If the chemical po- 
tential is tuned to /i = £^sb/2, the single particle band topol- 
ogy changes in a non-trivial way at this point (i.e. the first 
Chern number defined for two-band Bogoliubov-de Gennes 
Hamiltonian changes by 2) yielding an interesting phase di- 
agram which is similar to two coupled Majorana chains. In 
this case, however, instead of originating from physically dif- 
ferent chains, Majorana fermions emerge here from different 
transverse subbands. We first analyze the phase diagram as- 
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suming the fermion parity in each subband is preserved, and 
then discuss how various perturbations including interactions 
affect the topological phase diagram. 

We now provide a qualitative discussion of the topological 
phase diagram. For the moment, we neglect the band-mixing 

(12) (12) 

term Hso • As shown below, Hso drops out from the effec- 
tive low-energy theory at the multicritical point. Then, the two 
subbands decouple and one can define fermion parity in each 
subband, and one can introduce a Z2 topological invariant M 
(Majorana number)^^^for each subband Mi and M2 where 
M is defined as 



system to a low energy subspace £ <C E^h ' 



M = sgn 



Pf5(0)]sgn[Pf5 



±1. 



(8) 



The antisymmetric matrix 5 in Eq.([8]) represents the Hamil- 
tonian of the system in the Majorana basis^^ For the spe- 
cial momenta P = 0, 7r/a the antisymmetric matrix B can be 
constructed by the virtue of particle-hole symmetry: B{P) = 
HBdG{P)U EH and |25l where unitary matrix U defines 
particle-hole symmetry of BdG Hamiltonian 



ei/BdG(p)e- 



-HBdG{-p) 



(9) 



with Q = UK being the anti-unitary operator. Here K de- 
notes complex conjugation. The change of the topological 
invariant signals the phase transition with the phase boundary 
given by \V,\ = V/^+Ag and \V,\ = ^(/i-^sbj^+Ag for 
Uy = 1 and Uy = 2 subbands, respectively. At a special point 
in the phase diagram /i = E^h / 2 and Vx = ^/E^^/A + Aq, 
two topological phases can coexist and the symmetry group 
is Z2 ^ Z2, see Fig.[T]:. Around this point there are four dis- 
tinct phases: non-topological (no Majorana modes), topologi- 
cal with Majorana fermions originating either from n^y = 1 or 
Uy = 2 subbands, and the last one with two Majorana modes 
localized on each end, see Fig. [T]:. Thus, multiband semi- 
conductor nanowires are interesting both from a fundamental 
and a practical point of view as they offer a possibility to in- 
vestigate interaction between various topological phases in a 
realistic experimental system. 

From now on we focus on the multicritical point, and 
present a simple explanation of the topological phase transi- 
tion by deriving an effective long wavelength model around 
it. The topological phase transition requires vanishing of the 
excitation gap in the system for the reconstruction of the en- 
ergy spectrum to occui^. The quasiparticle excitation gap at 
the phase boundary vanishes as E{p) ~ \p\ |5|, and, thus, 
one can understand the phase diagram by deriving an effec- 
tive model in the spirit ofk-p perturbation theory. The phase 
transition between various topological phases can be captured 
by studying the Dirac-like Hamiltonian with the mass term(s) 
M changing sign across the phase boundary, see Fig.[T]:. We 
first calculate exact eigenstates around the multicritical point, 
and then evaluate small corrections due to the deviation of 
the physical parameters away from this point. Assuming that 
these terms are small compared to /i, Aq, ~ ^sb. we per- 
form the following canonical transformation and project the 



(10) 

Here 7^/^ are right/left-moving Majorana operators originat- 
ing from the first (a = 1) and second (a = 2) subbands, respec- 
tively. In the transformation above we kept only low energy 
degrees of freedom and neglected high-energy modes. (Note, 
however, that these high-energy modes are necessary to satisfy 
canonical anticommutation relations.) The amplitudes u± are 
given by 



U± 



2V2 



After some algebra, one arrives at the following effective 
Hamiltonian valid in the vicinity of the multicritical point: 



a=l,2 



dx 



(11) 



where a = aAo / \/E^- 



4An and the mass terms Ma can 



be written in terms of the deviations from the multicriti- 
cal point SVx 

^sb/2)/yE^ 



Vx - ^ElJA + Ag and Sfi = E,^{fi 



•4Ag: 



Ml/2 = -sVx±Sj:i. 



(12) 



The topological phase transitions can be classified in terms 
of the sign change of the mass terms M^, see Fig. [T]:. It is 
clear that the phase with Mi ^2 > is trivial since it is adia- 
batically connected to Vx = limit. A sign change of one 
of the two mass terms corresponds to a topological transition 
to a phase with a single Majorana mode. Finally, when both 
mass terms change sign we have two Majorana modes per end. 
However, the latter state is unstable against perturbations that 
break fermion parity in the individual chain and ultimately lift 
the degeneracy by hybridizing the Majorana modes. 



III. EFFECT OF SYMMETRY-BREAKING TERMS. 



Let us consider now effect of the symmetry -breaking terms, 
which couple Majorana modes 7*^^^ and 7*^^^ on the topo- 
logical phase diagram. The presence of such terms reduces 
the Z2 (8) Z2 symmetry of the model introduced in Eq. ( pT| ) 
to Z2 where the two phases corresponds to a different total 
fermion parity of the system^^ We now consider all possible 
^-dependent (independent of x) fermion bilinear perturbations 
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respecting the symmetry ([3]): 



(13) 



Hi^^^ = l dx[Ai2{d^c^-d^c^)^A[^{c^d^-c^d^)^h.c.] 

(14) 



In the gauge when Aq is real, A 12 should also be real which 
is assumed below whereas the interband tunneling matrix el- 
ements ti and t2 can be in general complex, i.e. ti = \ti\e^\ 
Using Eqs. ([T0|, one finds that 



{dtCi - d^c^) -- 



Ao 



y4AgTE^ 



/ (1) (2) , (1) (2)\ .-1^. 

[Tl Tl ^Tr Tr ) (15) 

/ (1) (2) (1) (2)\ 

[Tr'Tl -Tl'Ir) (16) 



sb 



( (1) (2) 

[1l Tl 



( (1) (2) (1) (2) 



(1) (2) 



(1) (2) 

■Tl Tr 



(17) 



(18) 



The bottom term in the above equation is allowed by symme- 
try ([3]) but invokes triplet pairing which is not generated in the 
setup considered here. Substituting these results into Eq.(p3]), 



we arrive at 



^^(12) 



2i|ti|sinaiAo / (i) (2) , (i) (2)\ 

[il rL+rRrR) 

) 



dx , ^ ^ 

2i\t2\ cos a2Ao 



/ (1) (2) (1) (2)^ 

[Tl Tr -Tr Tl 



(19) 



) (2) 

Tl 



(1) (2)1 



(20) 



Thus, at the level of the effective low-energy theory, there are 
only two types of allowed symmetry-breaking terms which 
couple Majorana fermions in different subbands, and mod- 
ify the phase diagram at the qualitative level. The physi- 
cal origin of these terms requires breaking translational sym- 
metry along 7/-axis. (Note that the spin-orbit band-mixing 

(12) 

Hamiltonian Hso does not contribute to the effective long- 

(12) 

wavelength model since H^o corresponds to ^2 -perturbation 
with a2 = 7r/2.) In an experimental system, the semicon- 
ductor/superconductor interface can break translational invari- 
ance along ?/-axis. Alternatively, one may apply ^/-dependent 
electrostatic potential to generate band-mixing terms. 

We can now analyze the full Hamiltonian Hq^H 12, where 
H12 reads 



i^i2 = i\i j 



dx 

L 
L 

1X2 I dx 



(1) (2) 



(1) (2) 



(1) (2) , (1) (2)' 



(21) 
(22) 



and understand effect of the symmetry-breaking perturba- 
tions. The modified phase boundary can be obtained by look- 
ing at closing of excitation gap. The lowest branch of the 
excitation spectrum is given by 



E{p) 



(23) 



and vanishing of the excitation gap at Px 
phase boundary: 



defines a new 



(24) 



Thus, at J/i = there is a window —\/\\ + A2 < 5Vx < 
Y^A^ + A2 where topological phases from different subbands 
coexist corroborating the results of Ref . 24 We note that this 
scenario occurs as long as the Hamiltonian ([2]) respects the 
symmetry ([3]). The alternative scenario where two topolog- 
ically trivial phases hybridize is forbidden by the symmetry 
([3]). As emphasized in Ref. |24|, this effect is particularly im- 
portant for experimental realization of the topological super- 
conducting phase in semiconductor nanowires where chemi- 



r 



cal potential fluctuations pose serious constraint. Indeed, the 
topological phase around this region is to a large extent robust 
against chemical potential fluctuations which now have to be 
dji Esh to cause the transition into the nontopological state. 

One can notice that at Sfi = and Ai = the quasipar- 
ticle excitation gap also closes at a finite momentum \px \ = 



\2 — 6V^/a. However, closing of the gap in this case does 
not correspond to a topological transition which should be ac- 
companied by a change of the topological invariant. As fol- 
lows from Eq. ([8]), the topological invariant can change only 
at the particle-hole symmetric points Px = O^n/a. Therefore, 
only gap closing at these point corresponds to a topological 
phase transition. 
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bosonized Hamiltonian H becomes 



0.5 514 



FIG. 2. (Color online) a) Topological phase diagram for two band 
semiconductor model in the presence interband superconducting 
pairing (i.e. yi — A 12 and y2 — 0) which breaks fermion parity 
in the individual subband and leads to the hybridization between the 
Majorana modes originating from riy — l and ny — 2 subbands. As 
a result, the enhanced Z2 Z2 symmetry at the multicritical point is 
broken down to Z2, compare with Fig[T]:. b) Shift of the topological 
phase boundary caused by interactions. Repulsive/attractive interac- 
tion lead to the suppression/enhansement of the topological phase at 
the multicritical point. 



IV. INTERACTION EFFECTS. 



dx 



25/1 . ^ 2Ai 

H sm Zip cos 2(p- 

27ra 27ra 



2A2 



(26) 



where v, K slyq related to the microscopic parameters of the 
model: v = 



Here a is a cutoff in the problem, which is related to the 
momentum bandwidth A ~ 1/a |26|. To simplify Hamil- 
tonian ([26]) we shift 2(p ^ 2(p — tan~^(Ai/(5/i) and in- 
troduce dimensionless parameters yi = SVxa/v, and y2 = 
^/Xl-\-5jl'^a/v. Then, the Hamiltonian becomes ([26]) 



' 2lT 
2^2 



dx 

L 

sin 2(p- 



K{d,^f^K-\djf 



2^1 
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2A, 



(27) 



The effective Hamiltonian H^^ describes the physics around 
the multicritical point in the presence of interactions and vari- 
ous symmetry-breaking perturbations. In the rest of the paper, 
we analyze Eq. ( [27] ) and compute corrections to the topologi- 
cal phase boundary due to the presence of the interactions. 



We now study the effect of interactions on the topologi- 
cal superconducting phase. For simplicity, we consider short- 
range interactions given by the Hamiltonian 



Hint = U / dx:{c{cx^d{dx)::{clcx'^dldx')'- 

(25) 



where U = 



In the last line of Eq. ([25]) we have 



projected the Hamiltonian to the low energy subspace. Com- 
bining all the terms, the full Hamiltonian for the interacting 
system becomes H = Ho-\-Hi2 +i^int • One can see that at the 
multicritical point (assuming Ai = and A2 = 0), the Hamil- 
tonian H maps onto Thirring model which can be solved ex- 
actly using bosonization. Note that the interaction term ( [25] ) 
does not break fermion parity in the individual chain, and pre- 
serves Z2 (8) Z2 symmetry present at the multicritical point. 
The effect of deviations from the multicritical point as well 
as effect of symmetry -breaking terms can be included in the 
model perturbatively assuming that these terms are small to 
begin with. We proceed using standard bosonization^^ by first 
introducing left and right-moving Dirac fermions tl^R/L 

yiR/L 



i7^y^)/V^, and then rewriting them in terms of 



A. Phase diagram for A2 = 0. 

It is instructive to first analyze the phase diagram with 
A2 =0. In this case, the Hamiltonian ( [27] ) appears in other 
systems such as classical two-dimensional XY model in a 
magnetic field and weakly coupled Heisenberg chains, see, 
e.g. Refs. ,27 - 29. and one can develop some intuition based 
on these analogies. 

Let us investigate how interactions modify the topological 
phase diagram of the system and study effects of the compet- 
ing relevant operators on the non-trivial criticality. Assuming 
that yi^2 and small interactions |K — 1| <C 1, one can 
obtain flow equations using a perturbative real-space RG ap- 
proaclT" 



dyi(l) 

dl 
dy2{l) 

dl 
d\nK 

dl 



{2-K)yr{l), 
(2-K-^)y2(l). 



K-^^-Kyl 



(28) 
(29) 
(30) 



the bosonic fields (p and 9: ipR/L 



The 



where / = ln[a/ao] is the flow parameter with ao being the 
initial value of the cutoff. Above RG equations reflect duality 
in the model: (p ^ at K ^ and yi ^ y2- One can 
see that in the vicinity of = 1, where we have a non-trivial 
critical point, both mass terms are relevant and flow to strong 
coupling under RG, i.e. each perturbation acting separately 
would yield a massive field theory. However, given that yi 
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and 7/2 couple to dual field operators corresponding to charge- 
density- wave pairing and Cooper-pairing, respectively, and 
drive the system to different ground states, the interplay be- 
tween them gives rise to a second-order phase transition at 
the intermediate coupling. Indeed, this is what happens at the 
self-dual point K = l. Away from this exactly solvable point, 
some intuition can be obtained by invoking scaling arguments. 
Using the analogy with 2D classical theory, the scaling theory 
of our zero-temperature ID problem can be easily formulated. 
Since the critical theory should be invariant under the rescal- 
ing, the singular part of the energy density satisfies the follow- 
ing scaling relations: 



fs[yi^y2] 



(31) 



The scale / can be fixed by requiring e^'^ ^ ^^^%/2 = 1- Sub- 
stituting /* into the energy density, one finds 



fs[yi^y2] =y| ' fs 



y2 



2-K 
2-K-^ 



yiA 



(32) 



The phase transition in this system occurs when 



2-K 
2-K- 



^2 " " yi ^ and, thus, the new phase boundary is 
given by 



when 25Q{l)a ~ 1. To the lowest order in yi, y2 and 5q, the 
RG flow of Sq is given by: 



(34) 



Since all perturbations yi , y2 and Sq are relevant, they flow to 
strong coupling under RG. If yi{l) reaches strong coupling 
before SQ{l)a becomes order one, the phase boundary can 
be calculated as before and is given by Eq. ( [33] ). The con- 
dition ^i(r) ~ 1 defines a new length scale This sce- 
nario is self-consistent as long as 2SQ{l*)a <C 1 which trans- 
lates into the requirement that SQ{0)a <C yi{0) 2-^ . In the 
opposite regime, the ?/i-term does not reach strong coupling 
whereas the RG flow of y2-term is unaffected by the A2 terms, 
and, thus, the system flows to an ordered phase determined 
by y2 cos{2(p). If y2 is exactly zero, the system is massless 
which corresponds to an incommensurate phase in the context 
of C-IC transition. For K = 1 this transition corresponds to 
closing of the gap at a finite momentum which occurs when 
1^2 1 > in agreement with the arguments presented 

above. As discussed in Sec. Illl| closing of the gap at finite 



momentum is not associated with a change of the topological 
invariant and, thus, C-IC transition in the bosonic language 
does not correspond to a topological transition. 



(33) 



This is one of the main results of the paper: the presence of 
weak interactions does not destroy the multicritical point but 
leads to a non-trivial renormalization of the phase boundary 
which now depends on the Luttinger liquid parameter K. One 
can easily see that Eq. ( [33] ) is consistent with the noninteract- 
ing (K = 1) result discussed earlier. In the case of a repulsive 
interaction (i^T < 1), the topological region at the multicritical 
point shrinks, see Fig. [2}). This result can be understood as a 
competition between the repulsive interaction and proximity- 
induced superconductivity, and in this regard, our conclusions 
are consistent with those of Refs. [3QH32j On the contrary, at- 
tractive interaction (K > 1) stabilizes the topological phase 
by expanding the area occupied by the topological phase, see 
Fig.|2b. 



B. Phase diagram for A2 / 0. 

Let us now consider effect of the A2-term on the topologi- 
cal phase diagram. It is helpful to first make a transformation 
20 ^ 20 ^ 2Sqx where Sq = KX2/V. First of aU, such a 
transformation eliminates dxO in Eq.([27]) and leads to spatially 
oscillating yi term. When ^2 = the Hamiltonian ( [27] ) corre- 
sponds to well-known Pokrovsky-Talapov model^^ which ex- 
hibits commensurate-incommensurate transition. If ^ga 1, 
the sm{20 + 2Sqx) is quickly oscillating and averages out to 
zero which reflects competition of yi and A2 terms. As a re- 
sult of this competition, the RG flow for yi has to be cutoff 



V. CONCLUSIONS. 

To conclude, we have studied interacting topological su- 
perconducting phases using a realistic model corresponding 
to the semiconductor nanowire in the limit of multiband oc- 
cupancy in contact with an s-wave superconductor. In the 
vicinity of the multicritical point in the phase diagram, the 
model considered here is equivalent to the one describing two 
coupled Majorana chains, and we have characterized the ef- 
fect of interactions as well as other perturbations on the topo- 
logical phase diagram. We find that moderate interactions do 
not affect the phase diagram qualitatively but lead to nontriv- 
ial quantitative changes in the phase boundary. Our results 
characterize the stability of the topological phase around the 
multicritical point against interactions and have important im- 
plications for the experiments on Majorana nano wires. Our 
results for the topological phase diagram around multicritical 
points are generic and should be applicable to situations when 
there is a larger number of channels occupied in the nanowire 
(i.e. more than two) as considered recently in Refl34] 
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